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ABSTRACT 

Spectropolarimetric observations of the pre-main sequence early-G star HD 141943 
were obtained at three observing epochs (2007, 2009 and 2010). The observations 
were obtained using the 3.9-ni Anglo-Austrahan telescope with the UCLES echelle 
spectrograph and the SEMPOL spectropolarimeter visitor instrument. The brightness 
and surface magnetic field topologies (given in Paper I) were used to determine the 
star's surface differential rotation and reconstruct the coronal magnetic field of the 
star. 

The coronal magnetic field at the 3 epochs shows on the largest scales that the 
field structure is dominated by the dipole component with possible evidence for the tilt 
of the dipole axis shifting between observations. We find very high levels of differential 
rotation on HD 141943 (~8 times the solar value for the magnetic features and ~5 
times solar for the brightness features) similar to that evidenced by another young 
early-G star, HD 171488. These results indicate that a significant increase in the level 
of differential rotation occurs for young stars around a spectral type of early-G. Also 
we find for the 2010 observations that there is a large difference in the differential 
rotation measured from the brightness and magnetic features, similar to that seen 
on early-K stars, but with the difference being much larger. We find only tentative 
evidence for temporal evolution in the differential rotation of HD 141943. 

Key words: Stars : activity - imaging - magnetic fields - Stars : individual : HD 
141943 



1 INTRODUCTION 

One of the key drivers of the solar dynamo is differential 
rotation. In the Sun strong shears occur in the interface layer 
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between the differentially rotating outer convective zone and 
the inner radiative zone, which rotates as a solid body. It is 
these shears that help convert the large-scale solar poloidal 
field into a strong toroidal component. However, for young, 
rapidly-rotating, solar-type stars a fundamentally different 
dynamo may be in operation. 

Spectr opolarimetric observations of yo u ng solar-type 
stars (i.e. 'Donati & Collier Cameron' '1997*; ^Donati et al.| 
19991: iD onati 1999; Donati ct al. 2003a; Marsdcn ct al. 200^; 
Dunstone et al.l 120081 ; IJeffers fc Donatil 120081 ') have shown 
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that their reconstructed magnetic topologies have large re- 
gions of near-surface azimuthal field. These are interpreted 
as the toroidal component of the large-scale dynamo field in 
the stars. The presence of these regions near the stellar sur- 
face has led to the belief that the dynamo operating in such 
stars is in fact distributed throughout the stellar convective 
zone, rather than being restricted to the interface layer as 
in the solar case. 

Differential rotation is still thought to play a role in the 
generation of magnetic fields in these stars, however how 
the dynamos in these stars operates is not well understood. 
Most theoretical models are based on our knowledge of the 
solar dynamo. The models of lKitchatinov fc Riidigerl (|l999l ) 
predict that the level of differential rotation on an early-G 
star should be greater than that on a mid-K star and that 
the level of surface differential rotation should decrease for 
stars with shorter rotat ional periods. While the models of 
iKiiker fc Riidigeil (|2005l also show that the level of differen- 
tial rotation on a star should be dependent upon its effective 
temperature (with hotter stars having higher levels of differ- 
ential rotation) but only weakly dependent on its rotation 
rate. These are predictions that we can now observationally 
test. 

There are several techniques we can use to observa- 
tionally measure the level of differential rotation on a star. 
Direct starspot tracking from multiple Dopple r images of 
the s urface features (i.e. ICollier Cameron. Donati fc Seme! 

cross-correlating two iiidepen dent Doppler images 
(i.e. Donati fc Collier Cameronlll997^ ■ incorporating a dif- 
feren tial rotation law i n to the Doppler imaging process 
(i.e. [Ponati et ah I [2OO0I : iPetit. Donati fc Collier CameronI 



Table 1. Fundamental parameters of HD 141943 from Paper I. 



Parameter 



value 



2002!) ■ or by Fourier analysis of stellar line profiles (i.e. 



Reiners fc SchmittI [200^ . 120031 ). 



Combining differential rotation measurements of a num- 
ber of young so lar-type stars , foun d using the Doppler 
imaging method, iBarnes et al.l (|2005D found that the level 
of differential rotation increases with stellar temperature 
with early-G stars having higher levels of differential ro- 
tation than lower-mass stars, in agr eement with the find- 
ings o f lKitchatinov fc Riidigeil (Il999l ) and Kiikcr fc R iidi geg 
II2OO3). The results also showed only a weak (if any) cor- 
relation with s tellar rotation rate, agai n in agreement with 
the findings of iKiiker fc Riidigeil (|2005l ). However, more re- 
cent mea surement s o f the young early-G sta r HD 171488 
llMarsden et al. 200^ : Ijeffers fc Donatil l2008l : iJeffers et all 
|2010D have shown this star to have an ev en higher level ofdif- 
ferential rotation than that indicated bv lBarnes et al.1 (I2OO5I ) 
with differential rotation measurements up to 10 ti mes the 
curren t solar value. This supports the findings of iReinerj 
l|2006f ) which show high levels of differential rotation on a 
number of F stars. 

Work bv lDonati. Collier Cameron fc Peti^ (|2003bl ') has 
shown that for some early-K stars the level of differential 
rotation measured from surface brightness features is consis- 
tently lower than that measured from the magnetic features 
(using the same dataset). This has led them to surmise that 
this is a result of the magnetic and brightness features be- 
ing anchored at different depths within the stellar convective 
zone and that the convective zone has a radially varying dif- 
ferential rotation structure, unlike the Sun. The work also 
showed that these early-K stars evidence temporal varia- 
tion in their levels of differential rotation that the authors 



Age 
Mass 

Photospheric temperature 
Spot temperature 
Unspotted luminosity 
Stellar radius 
iisini 

Radial velocity (I'rad) 

Inclination angle (i) 

Equatorial rotation period (Pcq) 



~17 Myrs 

~1.3 Mq 
5850 ± 100 

~3950 K 
2.8 ± 0.1 
1.6 ± 0.15 R0 
35.0 ± 0.5 km s" 
~0.1 km s~^ 

70° ± 10° 
~2.182 days 



"assumed errors (see Paper I). 

attribute to a feed-back mechanism in the stellar dynamo 
periodically converting magnetic energy into kinetic energy 
and vice-versa. 

As mentioned there is currently only one early-G star 
for which there are differential rotation measures from 
both brightness and magnetic features ob t ained at multiple 
epochs, HD 171488 (iMarsden et al.] I2OO6I : iJeffers fc Donatil 
l2008l : IJeffers et all I2OIOI ). This work has shown very little 
difference in the level of differential rotation measured from 
the magnetic and brightness features, although the measure- 
ment errors are large r than that found on the early-K stars 
llDonati et al ]|2003d ). Additionally, there appears to be little 
evidence of temporal evolution in the level of differential ro- 
tation on HD 171488. This has been speculated to be caused 
by the thinner convective zone of HD 171488 compared to 
that of the early-K stars previously studied. 

In order to expand the number of young early-G stars 
studied with spectropol arimetry this paper a long with the 
first paper i n the s eries iMarsden et al.l (|2010l . Paper I) and 
IWaite et al.l l|2010l ). presents differential rotation measure- 
ments, and magnetic maps, of two young early-G pre-main 
sequence (PMS) stars. Paper I deals with the reconstruction 
of the brightness and magnetic topologies of the star HD 
141943. This paper (Paper II) presents the coronal magnetic 
field reconstructions. Ha activity and differential rotation 
measurements of H D 141943. The third paper in the series 
IWaite et al.ll2010h deals with observations of HD 106506. 

As detailed in Paper I, HD 141943 is a young, active 
and rapidly-rotating PMS star. The stellar parameters have 
been determined in Paper I through the Doppler imaging 
process and are reproduced here in Table [1] 

The reconstructed brightness images show that it has 
a weak polar spot and a significant amount of low-latitude 
spot features at all epochs. The magnetic reconstructions 
show that its has a predominately non-axisymmetric radial 
field while its azimuthal field is predominately axisymmetric 
with ring of azimuthal field seen at the pole, similar to that 
of other active stars. 



2 OBSERVATIONS AND DATA REDUCTION 

HD141943 was observed at 3 epochs on the 3.9-m Anglo- 
Austrahan telescope (AAT), in March/ April 2007, April 
2009 and March/April 2010. A fourth observation set was 
taken in May 2006 (see Paper I) but due to the limited 
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nature of this set and the fact that it was only spectro- 
scopic observations we have not included it in our analysis 
here (Note: a differential rotation measure was attempted 
on this dataset but it was too limited to provide a re- 
sult). The three datasets used were spectropolarimetric ob- 
servations in left- and right-hand circularly polarised light 
taken using the University College London Echelle Spec- 
trograph (UCLES) and the SEMPOL sp ectropolarimeter 
ijSemel. Donati fc Reeslll993l : lDonati et a n i2003a ) visitor in- 
strument. 

The data were reduced using the ESpRIT (Echelle Spec- 
tra Red uction: an Int e ractiv e Tool) optimal extraction rou- 
tines of iDonati et al.l (|l997! ). We then used the technique 
of Least-Squares Deconvolution (LSD, iDonati et al.l Il997l l 
to sum the over 2600 photospheric spectral lines in each 
echelle spectrum in order to create a single high signal-to- 
noise (S/N) profile for each observation. Further details of 
the observations, including an observing log, and the reduc- 
tion process are given in Paper L 



3 RESULTS 

Brightness and magnetic maps of HD 141943 were created 
at the three epochs through the inversion of the observed 
Stokes I (brightness) and Stokes V (magnetic) LSD profiles. 
These images and the method used to create them are given 
in Paper L As an extension to these results we have mea- 
sured the differential rotation on the surface of HD 141943 
using these maps. We have also extrapolated the radial mag- 
netic field map to reconstruct the coronal magnetic field of 
the star and we have looked at the Ha emission from the 
star to look for prominence activity. 

3.1 Surface differential rotation 

If a star is observed for a number of rotations, then the 
surface features will evolve under the infiuence of the star's 
surface differential rotation. The surface differential rotation 
on a star can be determined from both its spot/brightness 
features as well as from its magnetic features. 

In order to measure the differential rotation on a star a 
simplified solar-like differential rotation law is incorporated 
into the imaging process: 

Q.{e) = ficq - dnsine{ra.A d"^), (1) 

where Q.{9) is the rotation rate at latitude 9, f^eq is the equa- 
torial rotation rate and dO, is the rotational shear between 
the equator and the poles (the differential rotation). The 
surface differential rotation is then determined by treating 
both fieq and dO. as free parameters and determining the 
best fit to the data using the x^ -minimisation meth o d. Thi s 
met hod is further des cribed in IDonati et all l|2000l . l2003bl ') 
and iPetit et al.l (|2002l ) . 

To determine the surface differential rotation usually re- 
quires that a significant part of the stellar surface is observed 
on at least two o ccasions preferably several days apart (see 
IPetit et al.ll2002l '). As the rotational period of HD 141943 is 
close to 2.2 days (see Table [!} this means that in order to 
determine the level of surface differential rotation the ob- 
servations need to cover a time base of around a week or 



Table 2. Surface differential rotation measurements for HD 
141943, taking into account the errors in the stellar parameters 
given in Table [Tl 



Epoch 








(rad d-l) 


(rad d-l) 




Stokes I (brightness) 




2010.244 


2.86 ± 0.02 


0.24 ± 0.03 




Stokes V (magnetic) 




2007.257 


2.88 ± 0.02 


0.36 ± 0.09 


2010.244 


2.89 ± 0.05 


0.45 ± 0.08 



more. We attempted to determine the level of surface differ- 
ential rotation on HD 141943 from the brightness and mag- 
netic features in each dataset. Unfortunately, the differen- 
tial rotation could only be measured from the March/ April 
2007 magnetic data and the March/April 2010 magnetic and 
brightness data. The other datasets did not produce a min- 
imum in the x^-landscape (see below) so the level of differ- 
ential rotation could not be determined. 

For the March/April 2007 magnetic data (Stokes V) the 
imaging code was forced to converge to a fixed magnetic field 
value of 91.3 G (with the level of magnetic field determined 
through an iterative process) for various values of Q,cq and 
dO.. This produced the reduced x^-landscape shown in the 
upper image of Fig. [T] Fitting a paraboloid to this data gave 
f^cq = 2.882 ± 0.009 rad d"^ and dQ, = 0.347 ± 0.035 rad 
d"'^, with the errors being la errors. 

For the March/April 2010 magnetic data (Stokes V) 
the imaging code was forced to converge to a fixed magnetic 
field value of 70.9 G and produced the reduced x^-landscape 
shown in the middle image of Fig. [1] Fitting a paraboloid to 
the central region of the data gave Slcq = 2.886 ± 0.009 rad 
d"^ and dO. = 0.439 ± 0.025 rad d'^ 

The level of differential rotation (dil) in 2010 for the 
magnetic features is higher that that found in 2007, with the 
difference corresponding to a ~ 2a change. This is a marginal 
change so we feel that it is only tentative evidence of tem- 
poral evolution in the differential rotation of HD 141943. 

For the March/April 2010 brightness data (Stokes I) the 
imaging code was forced to converge to a spot filling factor of 
0.029 (2.9 per cent) and produced the reduced x^-landscape 
shown in the bottom image of Fig. [T] Fitting a paraboloid 
gave: ficq = 2.856 ± 0.002 rad d~^ and dQ. = 0.232 ± 0.007 
rad d-\ 

For 2010 the level of differential rotation from the 
brightness features is significantly different to that found 
from the magnetic features, corresponding to a ~8(t change. 
This is a definite detection of different differential rota- 
tion rates for the brightness and magnetic features on HD 
141943. Fig. [2] shows a plot of the la error ellipses for all 
three differential rotation measurements on the same graph. 

These measures (and errors) of differential rotation as- 
sume that the stellar parameters used in the imaging code 
are correct and as shown by Table [T] they themselves have 
some uncertainty. The main parameters that will affect the 
measure of differential rotation are; thet;sini, inclination an- 
gle and rotational period. Taking these errors into account 
gave the surface differential rotation measures shown in Ta- 
ble [2] and also as the overplotted ellipses in Fig. [1] 
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Figure 1. Surface differential rotation x ™™i™isation for 
HD141943, March/April 2007 Stokes V data (upper plot), 
March/April 2010 Stokes V data (middle plot) and March/April 
2010 Stokes I data (bottom plot). The coloured images show the 
reduced values obtained from the maximum-entropy imaging 
code for various values of Qeq and dO assuming all stellar param- 
eters are correct. Darker regions correspond to lower reduced 
values. The coloured images project to itScr on both axes for the 
upper and middle plots and ±5(T for the lower plot. The overplot- 
ted ellipses show the la errors (projected onto both axes) for the 
differential rotation measures when taking into account the errors 
in the stellar parameters. 
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Figure 2. Plot showing the Icr error ellipses for the surface dif- 
ferential rotation for HD 141943. The upper ellipse is the Stokes 
V data from 2010, the middle ellipse is the Stoves V data from 
2007 and the bottom ellipse is the Stokes I data from 2010. These 
error ellipses have not taken into account the errors in the stellar 
parameters (see text). 



The level of differential rotation for the magnetic fea- 
tures on IfD 141943 is approximately 8 times the solar dif- 
ferential rotation rate. While the level of differential rotation 
for the brightness features is approximately 5 times the solar 
value. 

ft should be noted that the magnetic data for April 2009 
gave an extremely shallow paraboloid around fieq = 2.84 ± 
0.02 rad d"^ and dO. = 0.0 ± 0.1 rad d"\ However, the 
limited number of overlapping profiles and the poor quality 
of the observations early in the observing run (see Paper f) 
means that we believe this measurement of surface differen- 
tial rotation to be highly suspect. 

3.2 Coronal field extrapolation 

The coronal magnetic fields are extrapolated from 
the surface magnetogra ms using the "Potential Field 
Source Surface" method (lAltschuler fc Newkirk. Jr.1 Il969l : 
Ivan Balleeooiien et al.lll998l). Since the method has b een 



Jardine. Collier Cameron fc Donatil (|2002al ) we 



described in 

provide only an outline here. Briefly, the condition that the 
field is potential requires that (V x B = 0) and so we can 
write the magnetic field B in terms of a scalar flux func- 
tion ^ such that B = — V*!/. The condition that the field is 
divergence- free then reduces to Laplace's equation V^^* — 
with a solution in spherical co-ordinates (r, 6, (p) 



EE 

l — l m— — l 



-(' + !)■ 



Pim{e)e 



im(f) 



(2) 



where the associated Legendre functions are denoted by Pim- 
The coefficients ajm and bim are determined by imposing the 
radial field at the surface from the Zeeman-Doppler maps 
and by assuming that at some height Rs above the surface 
(known as the source surface) the pressure of the hot coronal 
gas overcomes the ability of the magnetic field to confine it. 
Thus, at the source surface the field lines are opened up to 
become purely radial, and hence Be{Rs) ~ B^{Rs) = 0. 

We determine the plasma pressure at every point by 
calculating the path of the field line through that point and 
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solving for isothermal, hydrostatic equilibrium along that 

path. In this case, the gas pressure is p = poe^/-^ 
where m is the mean particle mass, fcs is Bolzmann's con- 
stant, T is the temperature and p = po is the gas pressure 
at the base of the field line. We note that the integral in 
this expression is performed along the path of the field line 
and that Qs — (g.B)/jB| is the component of gravity along 
the field. We note that the plasma pressure is set to zero 
at any point where the field line through that point experi- 
ences a plasma pressure greater than the magnetic pressure 
somewhere along its length. In this case, we assume that this 
field line should have been forced open by the pressure of the 
plasma. The gas pressure at the footpoint of the field line 
po is a free parameter of this model. Following I Jardine et al.l 
(|2002al lbl) we choose to scale po to the magnetic pressure at 
the base of the field line, such that po = KBq where K is a 
constant that is the same on every field line. By scaling K 
up or down we can scale the overall level of the coronal gas 
pressure and hence the density and emission measure. Our 
model therefore has two parameters: the radius at which 
the field lines are opened up (the source surface) and the 
constant K which determines the gas pressure po at the 
base of each field line. These two parameters determine the 
magnetic field structure and X-ray emission measure of the 
closed-field regions of the stellar corona. These are given in 
Fig. for the three epochs. 



3.3 Hq variation 

The Ha line in active solar-type stars is often used as an 
activity indica tor with the line being "filled- in" in more ac- 
tive stars (i.e. ISoderblom et al]|l993l ). In addition, promi- 
nence activity in the stellar chromosphere can also be 
ma pped on such stars that show emission in the Hq lin e 
(i.e. ICollier Cameron fc Robinsonll 19891 : iDonati et a"l]|2000l ) . 
We have analysed the Ha line of HD 141943 in March/ April 
2007 and March/April 2010 as the most complete datasets 
to look for possible variations that may be attributable to 
prominence activity in the stellar chromosphere. As shown 
in Fig. 21 HD 141943 is a more active star than the Sun with 
a variable level of activity in 2007, but it does not have the 
overt Ha emission of a young T Tauri star. 

In order to look for prominence activity on HD 141943 
we have divided each of the Ha profiles in the March/ April 
2007 dataset by the mean Ha profile from the dataset. We 
did the same for the March/April 2010 dataset, dividing 
by the average Ha profile from the 2010 dataset. Dynamic 
spectra of these are displayed in Fig. [5] (2007 on the left 
and 2010 on the right) with darker areas showing regions of 
lower activity and lighter areas higher activity. 

As can be seen the Ha profile of HD 141943 in 2007 is 
variable in its activity level with regions of higher activity 
located between phases ~0.1 to ~0.6 while the other phases 
have lower activity than average. There would appear to 
be some evolution in the velocity of these act ive/non- active 
regions with most of the motion restricted to within the ro- 
tational velocity of the star, although there is some enhance- 
ment/reduction of the Ha emission outside these velocities. 
Given the inclination angle of HD 141943, prominences lo- 
cated at the co-rotation radius of the star are likely to be 
seen as absorption features quickly crossing the vsini range 




6555 6560 6565 6570 

Wavelength (A) 



Figure 4. Plot of the Ha profile of HD141943 at two epochs in 
March/ April 2007 along with the solar Ha profile taken with the 
same instrumental setup and shifted and convolved (by thejisini 
of HD 141943) to match the HD 141943 observations. 

of the star. The absorption features seen in Fig. [5] around 
phase '^O.SS appear to travel across the stellar profile, but 
not rapidly. A sine wave fitted to the peak/troughs of the 
Ha emission shows the amplitude of the sine wave to be re- 
stricted to within the vsini of HD 141943. Thus we believe 
that this region is located fairly close to the surface of the 
star and appears to be an inactive (in Ha) low-to-mid lati- 
tude feature on, or near, the stellar surface. Comparing this 
to the X-ray emission image (top-right image in Fig. [3]) we 
see that for phases around 0.0 there is little X-ray emission 
from the star. 

For 2010 there appears to be little or no variation in 
the Ha emission of HD 141943. The level of Ha emission at 
all phases in 2010 is similar to the highest level of emission 
in 2007. Thus it would appear that the Ha feature seen in 
2007 is not present on the star in 2010. 



4 DISCUSSION 

As mentioned in Paper I, HD 141943 is only th e second (or 
third including the results for HD 106506 by IWaite et al.l 
l2010h young early-G star for which the large-scale magnetic 
topology has been determined, the other being HD 171488 
(Marsdcn et al. 2006; Jcffcrs & Donati 2008; Jeffers et al.1 
l20ldh . In total there have been five young early-G stars for 
which differential rotation measures have been determined. 
The five stars are HD 141943, HD 1714 88, HD 106506, R5 8 
(|Marsden et al.ll2005al lbl) and LQ Lup (|Donati et al.ll200(]| '). 
The stellar parameters for all five stars have been given in 
Paper I. However, for ease of comparison we have given these 
in Table [3] and included the differential rotation results of 
the five stars. 



4.1 Surface differential rotation 

The differential rotation rate found for HD 14143 is one of 
the largest yet found using the Doppler imaging method and 
is similar to that of the other young early-G st ar HD 171488 
(|Marsden et al.1 l2006l : IJeffers fc Donatil l2008l : IJeffers et al] 
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Figure 3. The images on the left are the coronal magnetic structure of HD 141943 in 2007 (upper-left), 2009 (middle-left) and 2010 
(bottom-left) for a source surface of 4.8 R*. Closed field lines are shown in white while open field lines are shown in blue. The images 
on the right are the X-ray emission of HD 141943 in 2007 (upper-right), 2009 (middle-right) and 2010 (lower-right) for the same source 
surface and a coronal temperature of 2x10^^ K. All images are looking at phase 0.75 on the star. 



HD 141943 Halpha 2007 
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Velocity (km/s) Velocity (km/s) 

Figure 5. Dynamic spectrum of the variation in the Ha profile of HD141943, March/April 2007 (left) and March/April 2010 (right). 
All Ha profiles in each dataset have been divided by the average Ha profile of each dataset. The vertical thin lines represent the radial 
and rotational velocity of the star. 



Table 3. Comparison of the stellar parameters of the five young early-G stars that have had their surface differential rotation mea- 
sured using Dopple r imag ing. Except where noted, the data for HD 141943 come from this work and Paper I, that for HD 10 6506 
is from IWaite et al.l ll2010h and that for HD 171 488 fromlStra ssmeier et al.l l l2003f ), iMarsden et al.l ||2006|'). Ijeffers fc Donatl ||2008| 1 and 
Ijeffers et al.ll|2010l) . The data for R58 come fro mlM arsden e^a l. ( 2005a b|) and the data for LQ Lup from ^onat^^ al. (2000). The values 
for the depth of the convective zone are from lSiess et al. IHoOQ). For those stars with multiple measurements of differential rotation a 
value has been used with error bars large enough to encompass the entire range of the observations. 



Parameter 


HD 141943 


HD 106506 


HD 171488 


R58 


LQ Lup 


(B-V) 


0.65" 


0.605 


0.62" 


0.61* 


0.69'= 


Age (Myrs) 


~17 


<10 


30 - 50 


35 ± 5 


25 ± 10 


Mass (Mq) 


~1.3 


1.5 ± 0.1 


1.20 ± 0.02 


1.15 ± 0.05 


1.16 ± 0.04 


Radius (R0) 


1.6 ± 0.1 


2.15 ± 0.26 


1.15 ± 0.08 


1 10+0.17 
i-i»_o.io 


1.22 ± 0.12 


Inclination (°) 


70 ± 10 


65 ± 5 


60 ± 10 


60 ± 10 


35 ± 5 


Convective zone (R*) 


~0.16 [0.26 Rq] 


--0.22 [0.47 Rq] 


~0.21 [0.24 R0] 


~0.21 [0.25 R©] 


~0.26 [0.32 R0] 


Heq (rad d~^) Stokes I 


2.86 ± 0.02 


4.54 ± 0.01 


4.84 ± 0.14<* 


11.16 ± 0.04'' 


20.28 ± 0.01 


Qeq (rad d'^) Stokes V 


2.89 ± 0.05"^ 


4.51 ± 0.01 


4.81 ± O.IO'' 






dn (rad d"!) Stokes I 


0.24 ± 0.03 


0.21 ± 0.03 


0.41 ± 0.16^* 


0.08 ± 0.07'' 


0.12 ± 0.02 


dO, (rad d'^) Stokes V 


0.40 ± 0.13^* 


0.24 ± 0.03 


0.45 ± 0.06'' 







"from ICutispoto et aL I | |2002D : ''from lRandTciil 1I2OOII) . dereddened value; '^from lWichmann et al.l {l99^; ''multiple observations. 

It is also similar in level to tha t of the more ma- ferential rotation. However, the level of differential rotation 
ture planet-hostin g; late-F star Tau Boo l|Donati et al.ll200 for HD 141943 and HD 171488 are significantly above that 

iFares et al.ll2009l ) and is in agreement with the findings of of the other young early-G stars studied using the Doppler 

high lev els of differential rotation on inactive F stars by imaging method (see Table [31 . There appears to be little 

iReinerd I2OO6 ) using the line-profile method to measure dif- differences between the stars to explain the differences in 
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0.6 




-0.1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 
Convective zone depth (R*) 

Figure 6. Differential rotation {dD,) versus convective zone depth 
(as a function of st ellar radius) for young solar-type stars from 
I Barnes et al.l ll2005h and the new stars in Table [3] The c onvective 
zone depth has been determined from Sicss e t al.l l|2000h with an 
assumed error of ±0.05 R*. Dots show differential rotation mea- 
surements using brightness features and stars show differential 
rotation measures from magnetic features. For stars with mul- 
tiple measurements a value has been used with error bars large 
enough to encompass the entire range of the observations. 



differential rotation with the possible exception of their con- 

vect ive zone dept h. 

iBarnes et al.l (|2005h have shown that generally differ- 
ential rotation increases with stellar effective temperature, 
however, their dataset does not have any stars with a differ- 
ential rotation gre ater than dO ^ 0.2 r ad d~^. As an exten- 
sion of the work of lBarnes et al.l (|2005l ) we have plotted the 
differential rotation measured using the (Zeeman) Doppler 
imaging technique for these stars plus the additional stars 
listed in Table E] (HD 141943, HD 106506 and HD 171488). 
Fig. [S] plots this differential rotati on against stellar convec- 
tive z one depth (determined from lSiess. Dufour fc Forestinil 
I2OOOI ). 

As can be seen in Fig. |6]the level of differential rotation 
does appear to increase slightly with decreasing convection 
zones depth (with some scatter) until the convective zone 
depth reaches ^^0.2 (early-G stars) and then a dramatic 
increase in the level of surface differential rotation occurs. 
Why this should be is still not understood, but it appears 
that a change occurs in the ro tation of the conv e ctive zone 
at this depth. As suggested bv lJeffers fc Donatil (|2008l ). for 
stars with thinner convective zones, such as HD 141943, we 
may be seeing closer to the base of the convective zone and 
this could explain the higher levels of differential rotation 
seen on these stars^ 

Results from iDo nati et all l|2003bl ) have shown that for 
early-K stars the differential rotation measured from mag- 
netic features is higher than that measured from brightness 
features. They attribute this to the brightness and magnetic 
features being anchored at different depths in the stellar con- 
vective zone and the convective zone having a radi ally vary- 
ing di ffere ntial rotation (unlik e the Sun). Jeffers fc Donatil 
l|2008h and I Jeffers et al.l (|201G| ) found that there is virtually 
no change (within errors) in the differential rotation mea- 
sured from brightness and magnetic features for the early-G 
star HD 171488 although the errors in the differential rota- 



tion measurements are, for the most part, much larger than 
the errors for the early-K stars. Although their Stokes V dif- 
ferential rotation measurements are for the most part higher 
than those measured from Stokes L HD 106506 ( Wai te et al.l 
I2OIOI ) also shows only a small increase in the differential 
rotation from magnetic features over that from brightness 
features, but again the difference is within the errors of the 
measurements. 

In contrast to this our 2010 results for HD 141943 show 
a large difference between the differential rotation measured 
from the brightness features compared to that measured 
from the magnetic features, with the magnetic differential 
rotation being significantly higher than that from the bright- 
ness features. One possible reason for this difference could 
be the different latitude distributions of the brightness and 
magnetic features on HD 141943. In 2010 the spot features of 
HD 141943 are concentrated in a minor polar spot and some 
lower latitude features, while the magnetic features are more 
evenly distributed over the entire hemisphere (see Fig. 5 in 
Paper I). Such a difference in distribution may possibly lead 
to biases in the determination of the differential rotation, 
as the differential rotation is determined over a smaller (or 
larger) latitude range. Both HD 106506 and HD 171488 have 
larger polar spots and less low-latitude spot features than 
HD 141943. Thus, the two stars with more dominant polar 
spots would be expected to more affected by any bias in the 
latitude distribution between brightness and magnetic fea- 
tures, but neither star shows any significant disparity in dif- 
ferential rotation. The reason why HD 141943 alone shows a 
disparity in the differential rotation measurements between 
brightness and magnetic features remains unknown. A larger 
sample size is required to determine if HD 141943 is just a 
unique case. 

U nlike the results from early-K stars (iDonati et al.l 

l2003bl : I Jeffers et al.ll2007i ) the early-G star HD 171488 shows 
no evidence o f temporal evolution in its level of differen- 
tial rotation |Marsden et al ] I2OO6I : jjeffers fc Donatil I2OO8I : 
IJeffers et al.ll201Gl ) although again the errors are larger than 
the level of variation seen on early-K stars. Our magnetic 
results from 2007 and 2010 for HD 141943 show a slight dif- 
ference (see Section 13. ip but with only a ~2(j change this 
is only rather tentative evidence for temporal evolution in 
the star's differential rotation. The temporal variability seen 
on early-K stars is seen as evidence of a dynamo that pe- 
riodically converts magnetic energy into kinetic energy, and 
vise-versa, in the stellar convective zone. If such a mecha- 
nism is occurring in early-G stars (with thinner convective 
zones) it does not appear to have a large impact on the level 
of their differential rotation, although the errors in our dif- 
ferential rotation measures are too great to rule out a level 
of variability similar to that seen on early-K stars. 

Due to the rapid winding of magnetic field lines most 
likely to occur on stars with high levels of differential rota- 
tion, it is possible that such stars could have significantly 
shorter magnetic cycles than the 22-year solar magnetic cy- 
cle. Indeed the late-F star Tau Boo does in fac t appear to 
have a very short magnetic cycle (of ~2 years, iFares et al.l 
l2009h. but this sta r is also the host to a "Hot Jupiter" 
(jButler et al.l Il997l ) which may also be affecting its mag- 
netic cycle length. Both the early-G stars HD 141943 and 
HD 171488 have similar levels of differential rotation to that 
of Tau Boo and do not show any evidence of a magnetic po- 
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larity reversal over ^3 ye ars of observations (jjeffers et alj 
I2010I : iMarsden et al.|[2010l '). So the role of differential rota- 
tion in the magnetic cycles (or if they even have cycles) of 
these young stars is still unknown. 

4.2 Coronal magnetic field 

Fig. [3] shows the structure of the coronal magnetic field and 
also of the X-ray emission for the 2007, 2009 and 2010 mag- 
netic maps (see Paper I). Open field lines that would be 
X-ray dark and carrying the stellar wind, are shown blue, 
while the X-ray bright closed field lines are shown white. On 
the largest scales, the dipole component of the field domi- 
nates and the field structure for the three years is similar. 
The main change on the large scale is in the tilt of the dipole 
component of the field. Between 2007 and 2009, the dipole 
axis shifted from a latitude of 68 degrees to down to 52 de- 
grees and then back to 72 in 2010. On smaller scales, the 
field structure is very different between the three years and 
this is reflected in the X-ray images. As noted in Paper 1 
the poor quality of the 2009 dataset may play some role in 
these results. 

Fig. [7] shows the magnitude and rotational modula- 
tion of the emission measure, and also the emission-measure 
weighted density: 



/ nldV 
J nldV ' 



(3) 



These are calculated for a range of values of the base pres- 
sure, since without simultaneous X-ray observations we are 
unable to determine the value of base pressure. We also show 
in the left column results for a model with a small corona 
(extending to 2.5 R*) and in the right column results for 
a model with a larger corona (extending to 4.8 R*). It is 
immediately apparent that the small changes introduced by 
uncertainty in the size of the corona are small compared 
to the changes between years. We expect the corona to be 
populated with loops at a range of temperatures, but in 
the absence of any data that might determine the nature 
of the differential emission measure, we take the simplest 
approach. We therefore show results for a corona at a uni- 
form temperature of 2x10® K and for comparison, one at a 
uniform temperature of 2x10*^ K. 

As might be expected, at the lower temperatures, the 
densities and emission measures reach higher values. The 
range of values is however similar for all years and is typical 
for active stars. The largest difference between the three 
years is in the arrangement of the smaller scale loops and 
this results in a significant increase in the X-ray rotational 
modulation in 2009. 



5 CONCLUSIONS 

We have extrapolated the coronal magnetic field of HD 
141943 from surface magnetograms acquired at 3 epochs. 
These show that the large-scale field structure of the corona 
is dominated by a dipole component with the axis of the 
dipole shifting between the 3 epochs. The small scale struc- 
ture shows an increase in the modelled rotational X-ray 
modulation in 2009, compared to the other epochs (2007 

6 2010). 



The surface differential rotation of HD 141943 has been 
measured at 2 epochs (2007 & 2010) from both the surface 
brightness features (2010) and the surface magnetic features 
(2007 & 2010). The differential rotation measured from the 
magnetic features is one of the highest values of differential 
rotation measured on a young solar-type star and is simi- 
lar in level to the other young early-G star HD 171488, but 
higher than the more swollen star HD 106506. We thus con- 
clude that the depth of the stellar convective zone plays a 
strong role in the level of surface differential rotation seen 
on solar-type stars, with a large increase in differential ro- 
tation seen for star's with convective zone depths shallower 
than ~0.2 R*. 

The 2010 dataset for HD 141943 shows a large increase 
in the level of differential rotation measured from magnetic 
features to that measured from brightness features. This is 
similar to that seen on early-K stars but with a much greater 
difference and is in contrast to the results from other early- 
G stars which show little or no difference between the dif- 
ferential rotation measured from brightness and magnetic 
features. Our results only find tentative evidence for tem- 
poral evolution in the differential rotation of HD 141943. 
These results when combined with those from the early-G 
star HD 171488 (which shows no evidence of temporal evolu- 
tion in differential rotation) imply that early-G stars do not 
undergo large-scale evolution in their differential rotation. 
However, the errors in our measurements are too large to 
rule out small scale evolution in differential rotation similar 
to that seen on early-K stars. 

HD 141943 and stars of similar spectral type warrant 
further observations to determine what effect a shallow con- 
vective zone has on the differential rotation levels of such 
stars and indeed if they do show temporal evolution of their 
differential rotation as seen on early-K stars. 
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